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CN ■ Abstract 

^T*' We compute some components of the graviton propagator in loop quantum gravity, using the 

spinfoam formalism, up to some second order terms in the expansion parameter. 
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1 Introduction 



^ . An open problem in quantum gravity is to compute particle scattering amplitudes from the full 

■TJ" I background-independent theory, and recover low-energy physics p. The difficulty is that general 

'^ ■ covariance makes conventional n-point functions ill-defined in the absence of a background. A strategy 

^-p I for addressing this problem has been suggested in 0; the idea is to study the boundary amplitude, 

j^ i namely the functional integral over a finite spacetime region, seen as a function of the boundary 

\^^ ' value of the field (3] . In conventional quantum field theory, this boundary amplitude is well-defined 

f^ . (see ^ E] ) and codes the physical information of the theory; so does in quantum gravity, but in 

^^ I a fully background-independent manner |B]. A generally covariant definition of n-point functions 

^^. can then be based on the idea that the distance between physical points -arguments of the n-point 

J . ' function- is determined by the state of the gravitational field on the boundary of the spacetime 

bi)' region considered. This strategy was first implemented in the letter [7], where some components 

of the graviton propagator were computed to the first order in the expansion parameter A. For an 

implementation of these ideas in 3d, see |Hlinj. 

5^ ' Here we develop in more detail the calculation presented in of J2, and we extend it to terms of 

C^ I second order in A. We compute a term in the (connected) two-point function, starting from full non- 

perturbative quantum general relativity, in an appropriate large distance limit. Only a few components 
of the boundary states contribute to low order on A. This reduces the model to a 4d generalization 
of the "nutshell" 3d model studied in ^H]- The associated boundary amplitude can be read as the 
creation, interaction and annihilation of few "atoms of space" , in the sense in which Feynman diagrams 
in conventional quantum field theory expansion can be viewed as creation, interaction and annihilation 
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of particles. Using a natural gaussian form of the vacuum state, peaked on the intrinsic as well as 
the extrinsic geometry of the boundary, we derive an expression for a component of the graviton 
propagator. At large distance, this agrees with the conventional graviton propagator. 

Our main motivation is to show that a technique for computing particle scattering amplitudes in 
background-independent theories can be developed. (The viability of the notion of particle in a finite 
region is discussed in 11^. For the general relativistic formulation of quantum mechanics underlying 
this calculation, see ^21- On the relation between graviton propagator and 3-geometries transition 
amplitudes in the conventional perturbative expansion, see |13|.l 

We consider here riemaniann general relativity without matter. We use basic loop quantum gravity 
(LQG) results ^1 I15L [T^ , and define the dynamics by means of a spinfoam technique (for an intro- 
duction see ^lElEl and ^H^HI)- The specific model we use as example is the theory GFT/B, in 
the terminology of ^2], defined using group field theory methods j^|22]. On the definition of spin 
network states in group field theory formulation of spin foam models, see [23 and P21- The result 
extends immediately also to the theory GFT/C. These are background independent spinfoam theo- 
ries. The first was introduced in [21] and is favored by a number of arguments recently put forward 
1^1211 . The second was introduced in [2H] (see also |2Z1) and is characterized by particularly good 
finiteness properties |2H1- 

The physical correctness of these theories has been questioned because in the large distance limit 
their interaction vertex (lOj symbol, or Barrett-Crane vertex amplitude 1291 ) has been shown to include 
-beside the "good" term approximating the exponential of the Einstein-Hilbert action [301" also two 
"bad" terms: an exponential with opposite sign, giving the cosine of Regge action jSI] (analogous 
to the cosine in the Ponzano-Regge model) and a dominant term that depends on the existence of 
degenerate four-simplices |^ OSl 112] ■ We show here that only the "good" term contributes to the 
propagator. The others are suppressed by the rapidly oscillating phase in the vacuum state that peaks 
the state on its correct extrinsic geometry. Thus, the physical state selects the "forward" propagating 
Pl] component of the transition amplitude. This phenomenon was anticipated in jH5| . 

2 The strategy: two-point function from the boundary ampli- 
tude 

We begin by illustrating the quantities and some techniques that we are going to use in quantum 
gravity within a simple context. 

2.1 A single degree of freedom 

Consider the two-point function of a single harmonic oscillator with mass m and angular frequency 
uj. This is given by 

G^ihM) = {Q\x{h)x{t2m = {0\x e-i^(*i-*^) x|0) (1) 

where |0) is the vacuum state, x{t) is the Heisenberg position operator at time t and H the hamiltonian. 
We write a subscript o in 60(^1,^2) to remind us that this is an expectation value computed on the 
vacuum state. Later we will also consider similar expectation values computed on other classes of 
states, as for instance in 

G^{hM) = {'^M{tiXt2m- (2) 

Elementary creation and annihilation operator techniques give 

Go(ti,t2) = -^e-S-(*-*^). (3) 
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In the Schrodinger picture, the r.h.s. of Q) reads 



Go(ti,t2) = / dxidx2 ^po{xl) xi W[xi,X2]ti,t2\ X2 i/'o(a;2) (4) 



where ipoix) = (x|0) is the vacuum state and W[xi,X2]ti,t2] is the propagator, namely the matrix 
element of the evolution operator 

W[x^,X2-M,t2] - (a;i|e-*^(*i~*^V2>. (5) 

Recalling that 



^o(x)=^-e--^ (6) 

and (see for instance ^21) page 168 and errata) 

W{x,^X2;T)^J—^^e^^'-^^^^T (7) 

are two gaussian expressions, we obtain the two-point function (Q as the second momentum of a 
gaussian 

mUJ I 1 f ■ rn.u: {^f + ^l) <'°'"-T -2^i^2 ^^ 2 

Go{ti,t2) = —-\hr-. ^ dxidx2 X1X2 e' ^^ ='— ^ e " "^ , (8) 

Trn. V 2i sm uji J 

where the gaussian is the product of a "bulk" gaussian term and a "boundary" gaussian term. Using 

1 2tt 

dxidx2 X1X2 e"2(^^^) = A{1, (9) 

Vdet A 

the evaluation of the integral in ^ is straightforward. It gives 



„ . s rnijj / 1 27r , 1 , , 

G^{txM) ^ —r k\ -—. -, , A{} 10 



in terms of the inverse of the covariance matrix of the gaussian 






The matrix A is easy to invert and (|10() gives precisely Q . We will find precisely this structure of a 
similar matrix to invert at the end of the calculation of this paper. 

Notice that the two-point function (QJ can also be written as the (analytic continuation of the 
euclidean version of) the functional integral 

o{h,t2) = I Dx{t) x{h)x{t2) e^Ll^f--'*-/'^*). (12) 



Go 

where L is the harmonic oscillator lagrangian, and the measure is appropriately normalized. Let us 
break the (infinite number of) integration variables x{t) in various groups: those where t is less, equal 
or larger than, respectively, ti and ^2- Using this, and writing the integration variable x{ti) as xi and 
the integration variable ^(^2) as X2, we can rewrite H12|l as 



Go{ti,t2) = dxidx2 i^o{xi) xiW[xi,X2;ti,t2] X2 ■>jJo{x2) (13) 



where 



W[xi,X2-MM]^ I Dx{t)e^'2 ' ■ ' > (14) 



x{t2)=X2 



is the functional integral restricted to the open interval (ti,t2) integrated over the paths that start at 
a; 2 and end at a;i; while 



=(ti)= 



i?x(t)e'/-'~''^"''"/''^ 



Mx)= Dxit) e-l-'^ ' ' ' ' (15) 

Jx(-oc)=0 

is the functional integral restricted to the interval (— cx),ii). As well known, in the euclidean theory 
this gives the vacuum state. Thus, we recover again the form Q of the two-point function, with the 
additional information that the "bulk" propagator term can be viewed as the result of the functional 
integral in the interior of the (ii, ^2) interval, while the "boundary" term can be viewed as the result 
of the functional integral in the exterior. In this language the specification of the particular state |0) 
on which the expectation value of x{ti)x(t2) is computed, is coded in the boundary behavior of the 
functional integration variable at infinity: x{t) ^ for t — > ±cx). 

The normalization of the functional measure in (jl2|l is determined by 

i?x(t)e'/-"~''^^''"/'*'. (16) 

Breaking this functional integral in the same manner as the above one gives 



1 = / dxidx2 ^o(a;i) W[xi,X2;ti,t2] ipo{x2) (17) 

or equivalently 

l = (0|e-'^^(*i-*^)|0). (18) 

Let us comment on the interpretation of H13(l and p7(l , since analogues of these equation will play 
a major role below. Observe that 1)13(1 can be written in the form 

Ga{h,t2) ^ {Wt,,t2 |xi£2*o), (19) 

in terms of states and operators living in the Hilbert space )Cti.t2 = T^t ^ ^t2 (the tensor product of 
the space of states at time ti and the space of states at time ^2) formed by functions %p{xi,X2)- (See 
Section 5.1.4 of ^2] for details on Kt^^t^-) Using the relativistic formulation of quantum mechanics 
developed in [T^, thi s expression can be directly re-interpreted as follows, (i) The "boundary state" 
^a{xi,X2) = i}}{){xi)iI)q{x2) represents the joint boundary configuration of the system at the two 
times ti and ^2, if no excitation of the oscillator is present; it describes the joint outcome of a 
measurement at ti and a measurement at ^2, both of them detecting no excitations, (ii) The two 
operators xi and X2 create a ("incoming") excitation at i = t2 and a ("outgoing") excitation at t = ti; 
thus the state £10:2^0 can be interpreted as a boundary state representing the joint outcome of a 
measurement at ii and a measurement at ^2, both of them detecting a single excitation, (iii) The 
bra W4j^t2(j;i, 0:2) = VF[a;i,a;2;ii,i2] is the linear functional coding the dynamics, whose action on the 
two-excitation state associates it an amplitude, which can be compared with other similar amplitudes. 
For instance, observe that 

{Wt„t2 I i2 -^hm) = 0; (20) 

that is, the probability amplitude of measuring a single excitation at ^2 and no excitation at ii is zero. 
Finally, the normalization condition H17|l reads 

1 = {Wt,.t2 1*0); (21) 



which requires that the boundary state ^o is a sohition of the dynamics, in the sense that its projec- 
tion on ii is precisely the time evolution of its projection to ^2- As we shall see below, this condition 
generalizes to the case of interest for general relativity. We call (|21|l the "Wheeler-deWitt" (WdW) 
condition. This condition satisfied by the boundary state should not be confused with the normaliza- 
tion condition, 

l = (*o|*o), (22) 

which is also true, and which follows immediately from the fact that |0) is normalized in Tit. 

In general, given a state ^ e /Ct^^tj, the equations 

(m„t. I*) = 1; (23) 

and 

(*l*> = 1, (24) 

are equivalent to the full quantum dynamics, in the following sense. If the state is of the form 

\[/ = ■0f (g) -0;, then (123) and ^I^ imply that 

V-f = e-'^* Vi- (25) 

Finally, recall that a coherent (semiclassical) state i^i^ix) ~ er^i'^-i) +-kP^ is peaked on values 
q and p of position and momentum. In particular, the vacuum state of the harmonic oscillator is 
the coherent state peaked on the values g = and p — 0, with a = m to/h. Thus we can write 
'4'o = V'(g=o,p=o)- In the same manner, the boundary state ^q = ^o(a;i)V'o(a;2) can be viewed as a 
coherent boundary state, associated with the values gi = and pi = at ti and q2 — and p2 = at 
t2- We can write a generic coherent boundary state as 



'^qi,Pl.q2,P2{xi,X2) = ■4^(q^,p^){xi) "0(92 ,P2) (2^2 )• 



(26) 



A special case of these coherent boundary states is obtained when (gi,pi) are the classical evolution 
at time ii — ^2 of the initial conditions (92, P2)- That is, when in the ii — ^2 interval there exists a 
solution q{t),p{t) of the classical equations of motion precisely bounded by gi,pi, (72,^2, namely such 
that qi = q{ti),pi = p{ti) and (72 = (7(^2), P2 — p{i2)- If such a classical solution exists, we say that 
the quadruplet {qi,pi, 92, P2) is physical. As well known the harmonic oscillator dynamics gives in this 



case e-*^(*i-*^)*,,,p, = *g,,p,, or 



{Whm l*gi,Pi,«2,P2) = 1- (27) 

That is, it satisfies the WdW condition (|23|) . In this case, we denote the semiclassical boundary 
state a physical semiclassical boundary states. The vacuum boundary state ^0 is a particular case 
of this: it is the physical semiclassical boundary state determined by the classical solution q{t) = 
of the equations of motion, which is the one with minimal energy. Given a physical boundary state, 
we can consider a two-point function describing the propagation of a quantum excitation "over" the 
semiclassical trajectory q{t),p{t) as 

Gqi,p,,q2,P2(.tl,t2) = (V'(gi,pi)|a;(il)a:(i2)|V'fe,P2)> = (W^tl,t2 I ^1*2 *9i,pi,g2,P2)- (28) 

This quantity will pay a considerable role below. Indeed, the main idea here is to compute quantum- 
gravity n-point functions using states that describe the boundary value of the gravitatonal field on 
given boundary surfaces. 

There is an interesting phenomenon regarding the phases of the boundary state '^qi,pi,q2,p2 (2^1 j ^2) 
and of the propagator Wt-^^t2ixi,X2) that should be noticed. If pi and p2 are different from zero. 



they give rise to a phase factor e~'k'-P^^^~'P^^^\ in the boundary state. In turn, it is easy to see that 
M^ti,t2(^ii2;2) contains precisely the inverse of this same phase factor, when expanded around {qi,q2). 
In fact, the phase of the propagator is the classical Hamilton function Sti.t2{xi,X2) (the value of the 
action, as a function of the boundary values |12p. Expanding the Hamilton function around qi and 
g2 gives to first order 

dS dS 

Sti,t2{xi,X2) = S'ti.ta (91,92) + -^—{xi - qi) + ^—{x2 -92), (29) 

oxi 0x2 

but 

Giving a phase factor eH-(pi=^i-*P22:2)^ which is precisely the inverse of the one in the boundary state. In 
the Schrodinger representation of (|28|l . the gaussian factor in the boundary state peaks the integration 
around (91,(72); in this region, we have that the phase of the boundary state is determined by the 
classical value of the momentum, and is cancelled by a corresponding phase factor in the propagator 
W. In particular, the rapidly oscillating phase in the boundary state fails to suppress the integral 
precisely because it is compensated by a corresponding rapidly oscillating phase in W . This, of course, 
is nothing that the realization, in this language, of the well-known emergence of classical trajectories 
from the constructive coherence of the quantum amplitudes. This phenomenon, noted in |7] in the 
context of quantum gravity, plays a major role below. 

2.2 Field theory 

Let us now go over to field theory. The two-point function (or particle propagator) is defined by the 
(analytic continuation of the euclidean version of the) path integral {h = I from now on) 

Goix,y) = (O|0(a;)0(y)|O) = (O|0(f) e-»«(-o-^o) 0(y)|O) = J D^x) <P{x)<P{y) e'^I^l, (31) 

where the normalization of the measure is determined by 

1= f D(l){x) e^^I"^' (32) 

and the subscript reminds that these are expectation values of products of field operators in the 
particular state |0). These equations generalize (|12ll and (|15|l to field theory.^ As before, we can break 
the integration variables of the path integral in various groups. For instance, in the values of the field 
in the five spacetime-regions identified by t being less, equal or larger than, respectively, Xq and j/o- 
This gives a Schrodinger representation of the two-point function of the form 



Go{x,y) = / DipiDip2 ^poilpl) (pi{x) W[ipi,ip2;{xo - yo)] ^2(2/) ^^2)- (33) 

where ipi is the three-dimensional field at time ti, and ip2 is the three-dimensional field at time ^2- 
For a free field, the field propagator (or propagation kernel) 

I^((pi,<^2;T) = ((^i|e-*^^|^2). (34) 



'^A well-known source of confusion is of course given by the fact that in the case of a free particle the propagator JSJ 
coincides with the 2-point function of the free field theory. 



and the boundary vacuum state are gaussian expression in the boundary field ip — {ipi,ip2). These 
expressions, and the functional integral 1)33(1 . are explicitly computed in 5 . In a free theory, the 
boundary vacuum state can be written as a physical semiclassical state peaked on vanishing field and 
momentum tt, as in H26() : 



*o(¥'i, V2) = '^^i=o,Tn=o.,v2=o,7r2=o{'Pi,V2) = ■0o(<^i) ipo{'P2)- (35) 

Notice that the momentum n — -^ is the derivative of the classical field normal to E. 

More interesting for what follows, we can choose a compact finite region TZ in spacetime, bounded 
by a closed 3d surface S, such that the two points x and y lie on E. Then we can separate the 
integration variables in 1(31(1 into those inside 7?., those on E and those outside TZ, and thus write the 
two-point function ((31(1 in the form 

Go(x, y)= J D^ ip{x) ip{y) W[ip; S] ^olv?), (36) 

where if is the field on E, 

W[ip]Y.] ^ I D(j)R e-'^-^^''"'-^ (37) 

is the functional integral restricted to the region TZ, and integrated over the interior fields (jyn bounded 
by the given boundary field tp. The boundary state ^o('/') is given by the integral restricted to the 
outside region, TZ. The boundary conditions on the functional integration variable 

(^(x) -> 0, for |a;| -^ 00 (38) 

determine the vacuum state. In a free theory, this is still a gaussian expression in ip, but the covariance 
matrix is non-trivial and is determined by the shape of E. The state ^0 can nevertheless be still viewed 
as a semiclassical boundary state associated to the compact boundary, peaked on the value (^ = of 
the field and the value tt = of a (generalized) momentum (the derivative of the field normal to the 
surface) [Q. Equation 1(361) will be our main tool in the following. 

In analogy with 1(19(1 . equation l(3()() can be written in the form 

Go{x,y) = {W^\pix)ip{y)^o). (39) 

in terms of states and operators living in a boundary Hilbert space /Cs associated with the 3d surface 
E. In terms of the relativistic formulation of quantum mechanics developed in [12j . this expression 
can be interpreted as follows, (i) The "boundary state" ^0 represents the boundary configuration 
of a quantum field on a surface E, when no particles are present; it represents the joint outcome 
of measurements on the entire surface E, showing no presence of particles, (ii) The two operators 
(p{x) (p{y) create a ("incoming") particle at y and a ("outgoing") particle at x; so that the boundary 
state ip{x) (f{y)'i'o represents the joint outcome of measurements on E, detecting a ("incoming") 
particle at y and a ("outgoing") particle at x. (iii) Finally, the bra Ws is the linear functional coding 
the dynamics, whose action on the two-particle boundary state associates it an amplitude, which can 
be compared with other analogous amplitudes. The normalization condition for the measure, equation 
1(32(1 . becomes the WdW condition 

1 = {W^ 1*0), (40) 

which singles out the physical boundary states. 

Finally, as before, let q = {q,p) be a given couple of boundary values of the field ip and its 
generalized momentum on E. If there exists a classical solution (p of the equations of motion whose 



restriction to S is g and whose normal derivative to S is p, then we say that q = (QtP) are physical 
boundary data. Let \l/q be a boundary state in JCs peaked on these values: schematically 

*q((p)~e-/('^-«)'+*/^^. (41) 

If q = iq,p) are physical boundary data, we say that ^q is a physical semiclassical state. In this case, 
we can consider the two-point function 

G^{x,y)^{Wj:\0{x)0{y)^^) (42) 

describing the propagation of a quantum, from y to x, over the classical field configuration cj) giving the 
boundary data q == {q,p)- In the Schrodinger representation of this expression, there is a cancellation 
of the phase of the boundary state 'I'q with the phase of the propagation kernel Ws , analogous to the 
one we have seen in the case of a single degree of freedom. 

2.3 Quantum gravity 

Let us formally write H36|) for pure general relativity, ignoring for the moment problems such as the 
definition of the integration measure, or ultraviolet divergences. Given a surface S, we can choose a 
generalized temporal gauge in which the degrees of freedom of gravity are expressed by the 3-metric 7 
induced on S, with components jabix) a, 6 = 1, 2, 3. That is, if the surface is locally given by x^ = 0, 
we gauge fix the 4d gravitational metric field g^u(x) by 344 = l,(74o — 0, and ^ab — gab- Then the 
graviton two-point function H36|l reads in this gauge 

Gg'^'^lx, y) = J[Dj] h^'ix) h^%) W[t, n *o(^), (43) 

where h°-'^{x) = 7"''(x) — (5"^. As observed for instance in ISl, if we assume that W['y] S] is given by a 
functional integration on the bulk, as in H37|) . where measure and action are generally covariant, then 
we have immediately that WYl\ E] is independent from (smooth deformations of) S. Hence, at fixed 
topology of E (say, the surface of a 3-sphere), we have W['-^\ S] = W^[7], that is 

Gt'\x, y) = J[Dj] h'^\x) h^\y) W[^] ^o{^). (44) 

What is the interpretation of the boundary state ^0(7) in a general covariant theory? In the case 
of the harmonic oscillator, the vacuum state |0) is the state that minimizes the energy. In the case 
of a free theory on a background, in addition, it is the sole Poincare invariant state. In both cases 
the vacuum state can also be obtained from a functional integral by fixing the behavior of the fields 
at infinity. But in background-independent quantum gravity, there is no energy to minimize and 
no global Poincare invariance. Furthermore, there is no background metric with respect to which to 
demand the gravitational field to vanish at infinity. In fact, it is well known that the unicity and the 
very definition of the vacuum state is highly problematic in nonperturbative quantum gravity (see for 
instance 12,j), a phenomenon that begins to manifest itself already in QFT on a curved background. 
Thus, in quantum gravity there is a multiplicity of possible states that we can consider as boundary 
states, instead of a single preferred one. 

Linearized quantum gravity gives us a crucial hint, and provides us with a straightforward way to 
interpret semiclassical boundary states. Indeed, consider linearized quantum gravity, namely the well- 
defined theory of a noninteracting spin-2 graviton field /i^i, (x) on a flat spacetime with background 
metric 3°^,. This theory has a preferred vacuum state |0). Now, choose a boundary surface S and 
denote q — {q,p) its three-geometry, formed by the 3-metric Qab and extrinsic curvature field p"''', 



induced on I] by the flat background metric of spacetime. The vacuum state defines a gaussian 

boundary state on S, peaked around h = 0. We can schematicaUy write this state as ^s(^) ^ e J 
(In the conventional case in which E is formed by two parallel hyper-planes, the explicit form of this 
state is given in 221-) Now, on S there are two metrics: the metric q induced by the background 
spacetime metric, and the metric j = q + h, induced by the true physical metric 5^1, = g'^^, + h^i,, 
which is the sum of the background metric and the dynamical linearized gravitational field. Therefore 
the vacuum functional ^o(^) defines a functional ^q(7) of the physical metric 7 0/ E as follows 

*q(7)=*q(9 + /i) = *0(ft). (45) 

Schematically 

*q(7) = *o(/i) = *o(7 -q)^e- 1^''-'>^\ (46) 

A bit more precisely, as was pointed out in 7 , wc must also take into account a phase term, generated 
by the fact that the normal derivative of the induced metric does not vanish {q changes if we deform 
E). This gives, again very schematically 

*q(7)^e-/(^-«)'+*/p'' (47) 

as in 1)41(1 . Recall indeed that in general relativity the intrinsic and extrinsic geometry play the role 
of canonical variable and conjugate variable. As pointed out in |Zj, a semiclassical boundary state 
must be peaked on both quantities, as coherent states of the harmonic oscillator are equally peaked 
on q and p. The functional ^q of the metric can immediately be interpreted as a boundary state of 
quantum gravity, as determined by the linearized theory. Observe that it depends on the background 
geometry of S, because q and p do: the form of this state is determined by the location of E with 
respect to the background metric of space. Therefore (when seen as a function of the true metric 
7) there are different possible boundary states in the linearized theory, depending on where is the 
boundary surface. Equivalently, there are different boundary states depending on what is the mean 
boundary geometry g on E. 

Now, in full quantum gravity we must expect, accordingly, to have many possible distinct semiclas- 
sical boundary states ^q(7) that are peaked on distinct 3-geometries q = {q,p). In the background- 
independent theory they cannot be anymore interpreted as determined by the location of E with 
respect to the background (because there is no background!). But they can still be interpreted as 
determined by the mean boundary geometry q on E. Their interpretation is therefore immediate: they 
represent coherent semiclassical states of the boundary geometry. The multiplicity of the possible lo- 
cations of E with respect to the background geometry in the background-dependent theory, translates 
into a multiplicity of possible coherent boundary states in the background-independent formalism. 

In fact, this conclusion follows immediately from the core physical assumption of general relativity: 
the identification of the gravitational field with the spacetime metric. A coherent boundary state of the 
gravitational field is peaked, in particular, on a given classical value of the metric. In the background- 
dependent picture, this can be interpreted as information about the location of E in spacetime. In 
a background-independent picture, there is no location in spacetime: the geometrical properties of 
anything is solely determined by the local value of the gravitational field. In a background-independent 
theory, the dependence on a boundary geometry is not in the location of E with respect to a background 
geometry, but rather in the boundary state of the gravitation field on the surface E itself. 

Having understood this, it is clear that the two-point function of a background-independent theory 
can be defined as a function of the mean boundary geometry, instead of a function of the background 
metric. If q = {q,p) is a given geometry of a closed surface E with the topology of a 3-sphere, and 
\['q is a coherent state peaked on this geometry, consider the expression 

G^^'^^{x,y) - f[D^] h^\x) h^%) WVl] *q(7). (48) 



At first sight, this expression appears to be meaningless. The r.h.s. is completely independent from 
the location of S on the spacetime manifold. What is then the meaning of the 4d coordinates x 
and y in the l.h.s.? In fact, this is nothing than the usual well-known problem of the conventional 
definition of n-point functions in generally covariant theories: if action and measure are generally 
covariant, equation H31(l is independent from x and y (as long a,s x y^ y); because a diffeomorphism on 
the integration variable can change x and y, leaving all the rest invariant. We seem to have hit the 
usual stumbling block that makes n-point functions useless in generally covariant theories. 

In fact, we have not, because the very dependence of G^'^'^{x, y) on q provides the obvious solution 
to this problem: let us define a "generally covariant 2-point function" G'^'^'^{x,y) as follows. Given 
a three-manifold 5*3 with the topology of a 3-sphere, equipped with given fields q — {qab{^) , p'^'' {y)) , 
and given two points x and y on this metric manifold, we define 

G^^^'^(x, y) = J[Dj] /i"^(x) h^^iy) W^] ^^(j). (49) 

The difference between H48|) and (|49|l is that in the first expression x and y are coordinates in the 
background 4d differential manifold, while in the second x and y are points in the 3d metric manifold 
(5*3,5). It is clear that with this definition the dependence of the 2-point function on x and y 
is non trivial: metric relations between x and y are determined by q. In particular, a 3d active 
diffeomorphism on the integration variable g changes x and y, but also q, leaving the metric relations 
between x and y invariant. 

The physically interesting case is when q = {q,p) are a set oi physical boundary conditions. Since 
we are considering here pure general relativity without matter, this means that there exists a Ricci 
flat spacetime with 4d metric g and an imbedding "E : S3 —^ M , such that g induces the three metric 
q and the extrinsic curvature p on S3. In this case, the semiclassical boundary state ^q is a physical 
state. Measure and boundary states must be normalized in such a way that 



[Dj] W[j] v|/q(7) = 1. (50) 

Then the two point function ()49|l is a non-trivial and invariant function of the physical 4d distance 

L-dg(I](x),I](y)). (51) 

It is clear that if g is the fiat metric this function must reduce immediately to the conventional 2-point 
function of the linearized theory, in the appropriate large distance limit. 

This is the definition of a generally covariant two-point function proposed in |2], which we use 
here. 

Finally, the physical interpretation of (|49(l is transparent: it defines an amplitude associated to 
a joint set of measurements performed on a surface E bounding a finite spacetime region, where the 
measurements include: (i) the average geometry of E itself, namely the physical distance between 
detectors, the time lapse between measurements, and so on; as well as (ii) the detection of a ("outgo- 
ing") particle (a graviton) at x and the detection of a ("incoming") particle (a graviton) at y. The two 
kinds of measurements, that are considered of different nature in non-generally-relativistic physics, 
are on equal footing in general relativistic physics (see |T2j, pg. 152-153). In generally covariant quan- 
tum field theory, the single boundary state /i°^(a;)/i'^''(y)^q codes the two. Notice that the quantum 
geometry in the interior of the region TZ is free to fluctuate. In fact, W can be interpreted as the 
sum over all interior 4-geometries. What is determined is a boundary geometry as measured by the 
physical apparatus that surrounds a potential interaction region. 

Equation H49(l can be realized concretely in LQG by identifying (i) the boundary Hilbert space 
associated to E with the (separable jHSj) Hilbert space spanned by the (abstract) spin network states 
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s, namely the s-knot states; (ii) the hnearized gravitational field operators h°'^{x) and h'^'^{y) with the 
corresponding LQG operators; (iii) the boundary state ^q with a suitable spin network functional 
\l/q[s] peaked on the geometry q; and finally, (iv) the boundary functional VF[s], representing the 
functional integral on the interior geometries bounded by the boundary geometry s, with the W\s\ 
defined by a spin foam model. This, indeed, is given by a sum over interior spinfoams, interpreted as 
quantized geometries. This gives the expression 

Gf ^''(x, y) = ^ W[s] h'^\^) k\y) ^^[s]. (52) 

S 

which we analyze in detail in rest of the paper. The WdW condition reads 

1 = Y. ^w ^q[^]- (53) 

Using these two equations together, we can write 



a form that allows us to disregard the overall normalization of W and ^q. We analyze these ingredients 
in detail in the next section. 



3 Graviton propagator: definition and ingredients 

Equation (|52|l is well-defined if we choose a dynamical model giving W[s], a boundary state ^q[s] and 
a form for the operator h°-''{x). In the exploratory spirit of [2j, we make here some tentative choices 
for these ingredient. In particular, we choose the boundary functional W[s] defined by the group field 
theory GFT/B. We consider here only some lowest order terms in the expansion of W[s] in the GFT 
coupling constant A. Furthermore, we consider only the first order in a large distance expansion. Our 
aim is to recover the 2-point function of the linearized theory, namely the graviton propagator, in this 
limit. 

3.1 The boundary functional W[s] 

We recall the definition of W[s] in the context of the spinfoam theory GFT/B, referring to H^ and 
pn] for motivations and details. Wc follow the notation of (12j. The theory is defined for a field 
4> : (50(4))^ —> i? by an action of the form 

S[cf>] ^ S^iM + ^SinM. (55) 

The field <j) can be expanded in modes (jf'a^'/,''a^i (see eq. (9.71) of ,12 ). Notation is as follows. The 
indices j„, n = 1,...,4 label simple 5*0(4) irreducible representations. Recall that the irreducible 
representations of 5*0(4) are labelled by a pair of spins {j+,j-), corresponding to the split of so(4) = 
su{2) X su{2) into its self-dual and antiself-dual rotations; the simple representations are the ones for 
which j+ — j_ = J, and are therefore labelled by a single spin j. The index a„ labels the components 
of vectors in the representation j„ . The index i labels an orthonormal basis of intertwiners (invariant 
vectors) on the tensor product of the four representations jVi ■ We choose a basis in which one of the 
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basis elements is the Barrett-Crane intertwiner isci given in eq. (9.99) of |12j . Expanded in terms of 
these modes, the kinetic term of the action is (eq. (9.73) of [T^ l 



S'kin = 1^ Yl "^"l^ ^^Z 



CKti Jti,* 



The interaction term is (eq. (9.74) of |12| ) 



(56) 



(57) 



Here the notation is as foUows. The indices n and m run from 1 to 5, with n ^ m. Nnm = Nmn and 
^i^"* i = 4>'!i^''^^'']i''^^ i and so on. A(jnm,in) is the Barrett-Crane vertex amphtude. This is 



A{3n 



n^^ 



B{j. 



(58) 



where B{jnm) = -4(j„m,iBc) is the lOj symbol, given in eq. (9.102) of 521- In the following we use 
also the formal notation / 0^ = 5kin [0] and j (t>^ = Sint W\ ■ 



50(4)-invariant observables of the theory are computed as the expectation values 

1 ^" -/^^"^/^^ 



W\s\ = 



D0 fs{^) e~ 



(59) 



where the normalization Z is the functional integral without fs{4')^ and fs{<P) is the function of the 
field determined by the spin network s — (F, j/, i„). Recall that a spin network is a graph F formed by 
nodes n connected by links Z, colored with representations ji associated to the links and intertwiners 
in associated to the nodes. We note Inm a link connecting the nodes n and ttt,, and j„m = jmn the 
corresponding color. The spin network function is defined in terms of the modes introduced above by 



fs{<^) = llWz 



(60) 



Here n runs over the nodes and, for each n, the index m runs over the four nodes that bound the four 
links Inm joining at n. Notice that each index Unm = ctmn appears exactly twice in the sum, and are 
thus contracted. 

Fixed a spin network s, H59|l can be treated by a perturbative expansion in A, which leads to a 
sum over Feynman diagrams. Expanding both numerator and denominator, we have 



W[s] 



1 

Z'n 



D0/,(0)e-/*' 



(61) 



1 A 

1 A^ 

'Zo2{5iy 



D0/s(<^) / (b 



(/,»)e-/*'-l£*mi±:|0,/,We-K 



D0/s(<^) / ^ 



-U 



where Zq = jD(j>e J . As usual in QFT, the normalization Z gives rise to all vacuum-vacuum 
transition amplitudes, and it role is to eliminate disconnected diagrams. 

Recall that this Feynman sum can be expressed as a sum over all connected spinfoams a = 
{Ti,jf,ie) bounded by the spin network s. A spinfoam is a two-complex S, namely an ensemble of 



12 



Table 1: Terminology 



Od Id 2d 3d 4d 

Spin networks: node, link; 

Spinfoams: vertex, edge, face; 

Triangulation: point, segment, triangle, tetrahedron, four-simplex. 



faces / bounded by edges e, in turn bounded by vertices v, colored with representations jf associated 
to the faces and intertwiners ie associated to the edges. 

The boundary of a spinfoam a — (S,jj^,Ze) is a spin network s = (r,j;,i„), where the graph T is 
the boundary of the two-complex S, j; = // anytime the link I of the spin network bounds a face / of 
the spinfoam and i„ = ie anytime the node n of the spin network bounds an edge e of the spinfoam. 
See the Table ^ for a summary of the terminology. 

The amplitudes can be reconstructed from the following Feynman rules; the propagator 

'PiZ^'^,-k^' E n ^Or.A,,., '5„„„;^„, (62) 

ir(n) ". 

where 7r(n) are the permutations of the four numbers n = 1, 2, 3, 4; and the vertex amplitude 

^zr'" = ( n *»"»Bc) ( n '^"..."™.) ^o»™)' (63) 

n n^m 

where the index n = 1, ..., 5 labels the five legs of the five-valent vertex; while the index m ^ n labels 
the four indices on each leg. 

A Feynman graph has vertices v and propagators that we call "edges" and denote e. A spinfoams 
a is obtained from a Feynman graph by: (i) selecting one term in each sum over representations and 
one term in each sum over permutations (eq. (|62|l '). in the sum that gives the amplitude of the graph; 
(ii) contracting the closed sequences of (5a„a„ in the propagators, vertices and boundary spin-network 
function; and (iii) associating a face /, colored by the corresponding representation jf, to each such 
sequence of propagators and boundary links. See jl2| for more details. We obtain in this manner the 
amplitude 

^w = I E n ^i™(^/) n ^^(^n™) ( n<*"i*Bc) ) . (64) 

a,diJ=s fea- veer \nes / 

Here a are spinfoams with vertices v dual to a four-simplex, bounded the spin network s. f are the 
faces of a; the spins j^^ label the representations associated to the ten faces adjacent to the vertex v, 
where n ^ m = 1, ..., 5; dim(j) is the dimension of the representation j. The colors of a faces / of cr 
bounded by a link Z of s is restricted to match the color of the link: jf — ji. The expression is written 
for arbitrary boundary spin- network intertwiners i„: the scalar product is in the intertwiner space 
and derives from the fact that the vertex amplitude projects on the sole Barrett-Crane intertwiner. 
The relation between the different elements is summarized in Table |21 



13 



The sum (|64|l can be written as a power series in A 

oo 

W[s]=J2>^''Wk[s] 
with 



fc=0 



Wk\s] 






(65) 
(66) 



where a^ is a spinfoam with k vertices. 



Table 2: Relation between a triangulation and its dual, in the and 4d bulk and in its 3d boundary. 
In parenthesis: adjacent elements. In italic, the two-complex and the spin-network's graph. The 
spinfoam is cr = {/S.\,jf,ie)- The spin network is s = (A3,ji,z„). 



a: 



coloring 



Aa 



A? 



coloring 



4-siniplex 

tetrahedron 

triangle 

segment 

point 



vertex (5 edg, 10 fac) 
edge (4 faces) ie 

face jf 



tetrahedron node (4 links) i„ — ie 
triangle link ji = jf 

segment 
point 



Finally, recall that the last expression can be interpreted as the quantum gravity boundary am- 
plitude associated to the boundary state defined by the spin network s P2|- The individual spin 
foams a appearing in the sum can be interpreted as (discretized) spacetimes bounded by a 3-geometry 
determined by s. That is, (|64() can be interpreted as a concrete definition of the formal functional 
integral 

*[g] = / Dg e^^'^Rfg] 

Jdg=q 

where g is a 3-geometry and the integral of the exponent of the general relativity action is over the 
4-geometries g bounded by q. Indeed, lfi4|l can also be derived from a discretization of a suitable 
formulation of this functional integral. We now turn to the physical interpretation of this boundary 
3-geometry. 



(67) 



3.2 Relation with geometry 

In order to select a physically relevant boundary state ^q[s], we need a geometrical interpretation of 
the boundary spin networks s. To this aim, recall that the spinfoam model can be obtained from a 
discretization of general relativity on a triangulated spacetime. The discretization can be obtained as 
follows. 

We associate an R^ vector ef to each segment s of the triangulation. The relation with the 
gravitational field can be thought as follows. Introduce 4d coordinates x'^ and represent the gravita- 
tional field by means of the one-form tetrad field e^(x) — e^ {x)dx'^'' (related to Einstein's metric by 
9tj.v{x) = ej^{x)ei fj,{x)) . Assuming that the triangulation is fine enough for this field to be approxi- 
mately constant on a tetrahedron, with constant value e^, associate the 4d vector ef = e^Ax^ to the 
segment s, where Ax^ is the coordinate difference between the initial and final extremes of s. Next, to 
each triangle t of the triangulation, associate the bivector (that is, the object with two antisymmetric 
indices) 

(68) 



B{' ^ eiei 
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where s and s' are two sides of the triangle. (As far as orientation is kept consistent, the choice of 
the sides does not affect the definition of Bj.'^). Bj^"^ is a discretization of the Plebanski two-form 
j^ij — g/ ^ gJ _ ^jj^g quantum theory is then formaUy obtained by choosing the quantities Bj^"^ as 
basic variables, and identifying them with 5*0(4) generators J/"^ associated to each triangle of the 
triangulation, or, equivalently, to each face of the corresponding dual spinfoam. (For a compairaison 
with Regge calculus, see jSZ]-) 

The geometry is then easily reconstructed using the 50(4) Casimirs. In particular, the peculiar 
form H68|l implies immediately that 

eijKLB'/B^"^ ^ (69) 

any time t — t'oit and t' share an edge. Accordingly, the pseudo-scalar Casimir C — ^uklJI'^ Jt^^ = 
is required to vanish. This determines the restriction to the simple representations, which are 
precisely the ones for which C vanishes. 

The scalar Casimir C — ^Jl"' Jtij = \B^'^ Btij^ on the other hand, is easily recognized, using 
again l|^ . as the square of the area At of the triangle t. Indeed, calling ass' the angle between s and 
s', we have: 



1 



-{eiei - ef ef,)(es/es'j - e^jCs'i) 



2 * 2 

= Cs • Cs es' -Cs' - (cs • Cs'f = lesPl es'P (1- cos^ass') 
= (|es| \es'\ smass')^=Al (70) 

For simple representations, the value of C is j{j + 1). The quantization of the geometrical area, with 
j{j + 1) eigenvalues is of course a key result of LQG, reappearing here in the context of the spinfoam 
models. It is the LQG result that assures us that we can interpret it as a physical discretization and 
not an artifact of the triangulation of spacetime. 

An explanation about units is needed. Bj/^ has units of a length square, hence C has units [L]"^. 
In the quantum theory, B('' is identified with j/'^ and C has discrete eigenvalues. The identification 
requires evidently a scale to be fixed. This scale determines the Planck constant. A posteriori, we can 
simply reconstruct the correct scale by using again LQG, where the area eigenvalues are 



^. = ^^V7(7TT) (71) 

where 7 is the Immirzi parameter, which we fix to unit below, together with the speed of light c. This 
fixes the scale of the discretization (that is, it fixes the "size" of the compact 5*0(4) group in physical 
units). 

Next, consider two triangles sharing a side. Say the triangle t has two sides: the segments Si and 
S2 while the triangle t' has two sides Si and S3. Consider the action of the 5*0(4) generators on the 
tensor product of the representation spaces associated to the two (faces dual to the two) triangles. 
This is given by the operators j// = J/'^ + J^/ (we omit the tensor with the identity operator in 
the notation)). Equation 1)69(1 . for t ^ t' implies, with simple algebra, that the pseudo-scalar Casimir 
Ctt' — eijKbJu^ Ju'^ vanishes as well. This implies that the tensor product of the two representations 
associated with the triangles t and t' is -again- only allowed to contain simple representations. Let t 
and t' be two of the four triangles of a given tetrahedron. In the dual picture, they correspond to two 
faces joining along an edge e of the spinfoam. Then Ctt' is the pseudo-scalar Casimir of the virtual 
link that defines the intertwiner associated to this edge, under the pairing that pairs t and t' . The 
vanishing of Ctt' implies that this virtual link, as well, is labeled by a simple representation. In the 
model we are considering all internal edges are labeled by the Barrett-Crane intertwiner, whose key 
property is precisely that it is a linear combination of virtual links with simple representations for any 
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possible pairing of the four adjacent faces, thus consistently with Cu' = 0. This is in fact la raison 
d'etre of the Barrett-Crane intertwiner. 

Let us now consider the boundary s of the spinfoam a. A face / that cuts the boundary, labelled 
by a simple representation jf, defines a link / of the boundary spin network s, equally colored with a 
representation ji = jf. As we have seen, the quantity jf{jf + 1) is to be interpreted as the area of the 
triangle dual to the face /. This triangle lies on the boundary and is cut by the link I. 

Notice that we have precisely the LQG result that the area of a triangle is determined by the spin 
associated to the link of the spin network that cuts it. We can therefore identify in a natural way the 
boundary spin networks with the spin network states of canonical LQG. Recall that in LQG a basis of 
states of the quantum geometry of a 3d surface is labelled by abstract spin networks s. Since our aim 
here is not to fix the details of the physically correct quantum theory of gravity, but only to develop a 
general relativistic quantum formalism, we will do so in the following, disregarding some open issues 
raised by this identification (see below). 

The interpretation of the intertwiners at the boundaries is more delicate. Consider an edge e of 
a that cuts the boundary at a node n of s. The node n, or the edge e are dual to a tetrahedron 
sitting on a boundary. Let t and t' be two faces of this tetrahedron, and say, as above, that the 
triangle t has two sides Si and S2 while the triangle t' has two sides Si and S3. Consider now the scalar 
Casimir Cu' = J^^, Jw u on the tensor product of the representation spaces of the two triangles. 
Straightforward algebra shows that 

Cu' = \Ct\ + \Cv\ + 2nfnt.. (72) 

where nl — e^ jKLB^^t^ is and i^ is the normalized vector normal to t and t' (that is, to Si,S2 
and S3). Finally, fit ■ nf = At At' cos au' , where att' is the dihedral angle between t and t'. This 
provides the interpretation of the color of a virtual link in the intertwiner associated to the node, in 
the corresponding decomposition: if the virtual link of this intertwiner is simple, with spin jtt', we 
have 

jw Uw + I) ^ A'^t + ^t' + At At' cos att' ■ (73) 

That is, the color of the virtual link is a quantum number determining the dihedral angle cos att' 
between the triangles t and i'; or, in the dual picture, the angle between the two corresponding links 
that join at n. 

Once more, this result is exactly the same in 3d LQG. In this case, to each link is associated an 
SU{2) generator J*,i = 1,2,3, that can be identified with the SU{2) valued two-form _E* integrated 
on the dual triangle. The color of the link is the quantum number of the SU{2) Casimir C = 
(Jj + Jl,){Jti + Jfi)- Expanding, we have c = | Jtp + | Jfp + 2JtJt't or 

itt'iitt' + l)=A^t+Ar+ At At' cos att'. (74) 

where jtt' is the quantum number labelling the eigenspaces of C. We are therefore lead to identify 
the intertwiner ij , in the boundary spin network, with the intertwiner ij , in the LQG spinnetwork 
states, since they represent the same physical quantity. 

In fact, there is a key difference between H73|l and ()74|1 . In H73(l . jtf is the quantum muiiber labelling 
a simple 5*0(4) representation (recall 50(4) irreducibles are labelled by pairs of spins, which are equal 
for simple representations); while in l|74(l . ]tt' is the single spin labelling an SU{2) representation. Some 
potential difficulties raised by this difference are discussed in Appendix B. As argued in the Appendix, 
if we disregard these difficulties and we identify the intertwiner ij , with the LQG intertwiner ij , , 
we obtain simply and consistently 

{%„ Nbc) = i2jtt' + 1) = dirnUtt'). (75) 
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The details of this interpretation do not play a role in this paper. We leave a more complete discussion 
of this issue open. 

This completes the geometrical interpretation of all quantities appearing in the spinfoam model. 

3.3 Graviton operator 

The next ingredient we need is the graviton field operator. This is the fluctuation of the metric operator 
over the flat metric. At every point of the surface E we chose a local frame in which the surface is 
locally stationary: three coordinates x° with a = 1,2,3 coordinatize S locally, and the metric is in 
the "temporal" gauge: 1744 = 1,54a = 0. To the flrst relevant order, we define h°'^{x) = g°'^{x) — 5°'^. 
It is convenient to consider here the fiuctuation of the densitized metric operator 

h^\x) = {detg)g''''ix) - S"^ = E''\x)E'"{x) - S"''. (76) 

In the linear theory, the propagators of the two agree because of the trace- free condition. To determine 
its action, we can equally use the geometrical interpretation discussed above, or, directly, LQG. We 
study the action of this operator on a boundary spin network state: 

£;"(f)^"(x)|s). (77) 

Let us identify the point x with one of the nodes n of the boundary spin network s. Equivalently, with 
(the center of) one of the tetrahedra of the triangulation. Four links emerge from this vertex. Say 
these are e/, 1=1,2, 3, 4. They are dual to the faces of the corresponding tetrahedron. Let n^ be the 
oriented normal to this face, defined as the vector product of two sides. Then E{n)^^ = _E°*(x)ri^ can 
be identified with the action of the an SU{2) generator J* on the edge e/. We have then immediately 
that the diagonal terms define diagonal operators 

E''{n)E^{n)\s) = {S^fiGf jaija + l)\s) (78) 

where ja is the spin of the link in the direction a. The non-diagonal terms, that we do not consider 
in the following, are given in Appendix C. 

3.4 The boundary vacuum state 

As discussed in Section 2, the propagator will depend on a geometry q of the boundary surface 
E. Let us begin by choosing this 3d geometry. Let q be isomorphic to the intrinsic and extrinsic 
geometry of the boundary Eq of a 4d (metric) ball in Euclidean i?** with given radius, much larger 
than the Planck length. We want to construct the state ^li'q[s]. (On the vacuum states in LQG, see 
123 El 11011111112 123 •) Below we shall only need the value of ^I'q[s] for the spinnetworks s = (TJiJn) 
defined on graphs F which are dual to 3d triangulations A. We identify each such A with a fixed 
triangulation of Eq. 

We assume here for simplicity that, for each graph, ^q[s] is given by a function of the spins of s 
which is non- vanishing only on a single intertwiner on each node, which prjects on the ibc intertwiner 
under H75() . This will play no role in this paper, because, as we shall see, we compute only diagonal 
components of the propagator, which do not depend on the intertwiners. The precise role of the 
intertwiners, and other choices for the intertwiner dependence of the boundary state, will be discussed 
elsewhere. 

The area Ai of the triangle t; of A, dual to the link I, determines background values j/"' of the 
spins j;, via 



Ai^SttUG ^jrHj{°^ + l). (79) 
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We take these background values large with respect to the Planck length, and we will later consider 
only the dominant terms in l/j/°'. 

We want a state ^q[s] — \I>q(r, j), where j — {ji}, to be peaked on these background values. The 
simplest possibility is to choose a Gaussian peaked on these values, for every graph T 



*qH . Cr exp -± V au> '-^-^ ^^— V + ^ Y. '^f 3i (80) 




where I runs on links of s, an is a given numerical matrix, k G (0,2) (see below), and Cr is a 
graph-dependent normalization factor for the gaussian. 

The phase factors in (|80|l play an important role [7]- As we know from elementary quantuin 
mechanics, the phase of a semiclassical state determines where the state is peaked in the conjugate 
variables, here the variables conjugate to the spins j;. Recall the form of the Regge action for one 
simplex, S'Rcggc = X)/ ^i{ji)ji^ where $z(j;) are the dihedral angles at the triangles^, which are function 
of the areas themselves and recall that dS^cggc/ dji = $/. It is then easy to see that these dihedral 
angles are precisely the variables conjugate to the spins. Notice that they code the extrinsic geometry 
of the boundary surface, and in GR the extrinsic curvature is indeed the variable conjugate to the 
3-metric. Thus, $; are determined by the dihedral angles of the triangulation A. 

Concerning the quadratic term in H80|l we have put the (l/j/"')''^^ factors in evidence because we 
want a semiclassical state for which the relative uncertainties of area and angle become small when 
all the areas are large, namely in the large distance limit in which all the spins jV' are of the order of 
a large jl. That is, we demand that 

AA A$ 

— >0 and -j- -> 0, when j/"> - jl -^ oo. (81) 

Assuming that the matrix elements CK(i)(i') ^ a do not scale with j^, the fluctuations determined by 
the gaussian state (P^ are of the order 

A,^|^, A^^p (82) 



Therefore, since angles do not scale. 



AA A, ,1/-^ 



A$^^. (83) 



(|81|l and (|83|l restricts to fc G (0, 2). From now on, we choose fc = 1. That is 



*.[.] - Cr exp -± V au, ^^^^ -^^^ + ^ V ^J"^ j, ■ (84) 




The need for this dependence on the scale of the background of the covariance matrix of the 
vacuum state was been pointed out by one of us in the 3d context ^ and by John Baez in the 4d 



^ These are angles between the normals to the tetrahedra, and should not be confused with the angles between the 
normals to the faces, which are related to the intertwiners, as we discussed in Section 3.2. 
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case, following numerical investigation by Dan Christensen and Greg Egan, that have shown that in 
the absence of this dependence the width of the gaussian is not sufficient for the approximation taken 
above to hold ^^. 

A strong constraint on the graph-dependent constants Cr and matrix aw is given by the WdW 
condition (|53|l . which requires the state to satisfy the dynamics. The physical interpretation of the 
matrix aui is rather obvious: it reflects the vacuum correlations, and is the analog of the covariance 
matrix in the exponent of the vacuum functional in the conventional Schrodinger representation of 
quantum field theory. The physical interpretation of the Cr coefficients is less clear to us; it bears on 
the core of the diff-invariant physics of loop quantum gravity, and will be discussed elsewhere. 

3.5 The lOj symbol and its derivatives 

Baez, Christensen and Egan have performed in \IM\ a detailed numerical analysis, which has lead 
them to conjecture that if we rescale all spins by a factor A, then for large A the lOj symbol can be 
expressed as a sum of two terms. 



BU^J) = Yl -^('^) '=°^ \SRegge{<j) + kj] + D(j,fc). (85) 



4J 

a 

P{a) is a slowly varying factor, that grows as A~^/^ when scaling the spins by A. To understand this 
formula, consider a 4-simplex in i?^, with triangles ty having areas Aij = \Jjij{jij + !)• In general, 
there may be several distinct 4-simplices with triangles having these areas: let's label the distinct 
4-simplices with a discrete label a. Each triangle Uj separates two boundary tetrahedra Ti and r, of 
the 4-simplex. Each tetrahedron Ti defines a normal vector n^, normalized and normal to all sides 
of the tetrahedron. The angle $ij between the normals Ui and rij is the dihedral angle between 
the tetrahedra Ti and tj . (The triangles tij are in one-to-one correspondence with the links I of the 
boundary spin network, hence the notation $y is consistent with the notation $/ used above.) For a 
fixed a, we can compute the dihedral angles $ij as a function ^ijijij), of the areas Aij, hence of the 
10 spins jij. The Regge action associated to the 10 spins is 

5'Regge (ct) = J2 ^v ^ij iJij ) • (86) 



It is characterized by the fact that 



^%f^-<i>.fe); (87) 



that is, the derivative with respect to the jij in the angles does not contribute to the total derivative 
(this is the discrete analog to the fact that when we vary the Einstein-Hilbert action with respect to 
the metric, the metric variation of the Christoffel symbols does not contribute.) 

The form H85|l for the lOj symbol has later been confirmed by detailed analytical calculations by 
Barrett and Steele |S2I and by Freidel and Louapre [HSI- As first noted in 30 , the first term in 185|l 
is very good news for quantum gravity: it indicates that the lOj symbols are indeed related to 4d 
general relativity. 

On the other hand, to understand the origin of the second term D(jik) in H85|l . recall that the lOj 
symbols can be expressed in the form 

B{hj)= dy, YY . (88) 
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Here ^ij is the angle between the two unit vectors yi and yj . The large j behavior of this expression 
has been evaluated in |22| and |8Hj using a stationary phase approximation. It turns out that the 
integration variables yi admit a very interesting geometrical interpretation as the normals to the 
tetrahedra yi = rii. The stationary points of the integral can be interpreted as different geometrical 
configurations. Some stationary points are given by non degenerate tetrahedra. These yield the Regge 
term in l)85|l . But there are also contributions to the integral coming when two of the yi are parallel, or 
more in general when the linear span of the five yi is of dimension smaller than four. These degenerate 
contributions yield the D{jik) term. 

The bad news is that this degenerate term strongly dominates for large j. This fact casts a thick 
shadow of doubt over hope that a Barrett-Crane spinfoam model could yield the correct general 
relativistic dynamics. In fact, the discovery of the degenerate contributions is one of the sources 
of a recent decrease in interest in these models. However, light comes back into the shadow, in 
consideration of the results of the previous section. In fact, observe that what enters in the expression 
we have found in the previous section is not the lOj symbol itself, but rather a second derivative 
of the lOj symbol, because of the field insertions in (|54|l . The fact that the degenerate term D{jij) 
dominates over the Regge term at large j does not imply that its second derivative dominates as well. 

A hint that this hope may be correct comes from the following naive argument. Degenerate 
contributions arise when the denominator in the integrand of H88(l vanishes. The formal second 
derivative with respect to jij, considered as a continuous variable, is 

dja 7(53)5 fjr sm($y) 

which decreases the order of the divergence in $y =0. A better version of this argument, due to 
Greg Egan 0S]j is the following. The diagonal term of the discrete second derivative is i?(j + 1, ...) — 
2B{j, ...) + B{j — 1, ...). Using the trigonometric identity: 

sin((a + 2)t) - 2 sin(ai) + sin((a - 2)t) ^ -4 sm{t)^ sin(at) (90) 

we see that the effect of taking the diagonal discrete second derivative is to multiply the relevant 
kernel by — 4sin(<I>)^. This should then eliminate all the fully degenerate points from the integral. 
(The elimination of the fully degenerate points does not necessarily mean that the only remaining 
contribution is the Regge one, since, as shown in 32^, there is a complex zoology of other degenerate 
contributions to the integral. Numerical analysis to address some of this issues is being developed by 
Dan Christensen 1461.) 

4 Order zero 

Let us begin by evaluating the general covariant 2-point function to order zero in A. To this order 

Wo[s] = Z^^ lDcl,fs{^)e-I*\ (91) 

The Wick expansion of this integral gives non-vanishing contributions for all s with an even number 
of nodes. Since there are no vertices, each of these contributions is simply given by products of face 
contributions, namely products of dimensions of representations. The 2-point function H52() reads 

G^^^'^(x,y) - ^W^oW h'^\^)h^''{y)^^[.s]. (92) 
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Inserting (|84|l . we have 
Gf '^(x,y) =J2Wo[s] h^\^)h^'iy) Crexpj-ij^a,. ^^ Zi:_|^ + , ^ $(») , A . (93) 

s I W VJI VJl' I ) 

We are interested in this expression for large jY' ■ In this regime, the gaussian effectively restricts the 
sum over (a large region of) spins of order f°\ Over this region, the phase factor fluctuates widely, 
and suppresses the sum, unless it is compensated by a similar phase factor. But Wo[s] contains only 
powers of ji 's, and cannot provide this compensation. Hence we do not expect a contribution of zero 
order to the sum. The only exception can be the null spin network s = 0, which gives iy[0] = 1 
because of the normalization. Hence, to order zero 

G^^^'^(x,y) = Worn h'^'i^yh'-'iy) C0. (94) 

But is reasonable to assume that the semiclassical boundary state on a macroscopic geometry q has 
vanishing component on s = 0, whose interpretation is that of a quantum state without any volume. 
Hence we take C0 = 0, and we conclude that the 2-point function has no zero order component in A. 

This result has a compelling geometrical interpretation. The sum over spinfoams can be interpreted 
as a sum over 4-geometries. The boundary state vE'q[s] describes a boundary geometry which has a 
nontrivial extrinsic curvature, described by the phase of the state. In the large distance limit, we expect 
semiclassical histories to dominate the path integral. These must be close to a classical solution of the 
equations of motion, fitting the boundary data. Because of the extrinsic curvature of the boundary 
data, it is necessary that the internal geometry has non vanishing 4- volume. A round soccer ball must 
have volume inside. But the four-volume of a spinfoam is given by its vertices, which are dual to the 
four-simplices of the triangulation. Absence of vertices means absence of four-volume. It is therefore 
to be expected that the zero order contribution, which has no vertices, and therefore zero volume, is 
suppressed by the phases of the boundary state, representing the extrinsic curvature. 

Let us therefore go over to to the first order in A. 



5 First order: the 4d nutshell 

Consider a spinfoam a, dual to a single 4-simplex. Its the boundary spinnetwork has five nodes, 
connected by 10 links, forming the five-valent graph T^. That is 




(95) 



334 



The boundary function fs{(t>) determined by this spin network is 

f ((h) ^ \ ' J)'^12Q;i3Q;i4«i Jya2ia230i24,i2 Jya3ia32aaii3 iai2ai3aiiii J^a5ia52a5ai5 fQfil 

JsyHJ) / ^ 'rj-^2Jl3Jlijl5 ^J21J23J24J25 ^J31 J32j34 J35 ^J4lji2Ji3J4,5 ^J5lJ52J53J54 ^ ' 
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where a„m = a„i„. This can be compactly written as 

fM = E n Cr" (97) 

anjn n— 1,5 

This is a monomial of order five in the field, and is an observable in the group field theory. Its 
expectation value is given by H59|l . We consider the perturbative expansion (|62|l . At order A, the only 
term remaining is 

M/ W = A /■ D(/) fs {<t>) ( I A e-I'^\ (98) 



The Wick expansion of this integral gives one vertex v and five propagators 

A 
5! 



^'w --Ull ^^^.fc^O^fc:^;' (99) 

n=1.5 



where repeated representation indices are summed over. This expression still contains many terms due 
to the summation over the permutations in the propagators. Recall that we can give the geometrical 
meaning of a face to each closed sequence of deltas in this expression. Each face contributes with a 
factor equal to the dimension of the representation. The dominant term for large representations is 
therefore the one with the largest number of surfaces. A short reflection will convince the reader that 
this is the term in which the surfaces correspond precisely at the faces of the dual of a four-simplex. 
That is, the dominant term of W[s] to order A is 



^ (n(*"l*Bc)j ( n dim(^™")) ^(J-™)- 



^W - ^ M1(*«Nbc) 11 dim(j„,„) B{j:J. (100) 

\ n / \7i<?n / 

Since we have chosen boundary states peaked on an intertwiner that projects on zbci this reduces to 

(101) 



^w = M n dini(j,„„) ) ^(j-;;™)- 

\n<m / 



This is the dominant term of the connected component of the amplitude for the boundary spin network 
considered, in the limit of large representations. This is the expression we will use within equation 

The value of ^q[s] on the spin-networks s — (T^^jnm) (here n,m = 1,...,5) can be determined 
by triangulating Eq with the 3d triangulation formed by the boundary of a regular four-simplex 
of side L. The area of the triangles is Al — V5L^/4:. Then H79|l implies that jiim — Jl where 
87r?iG\/ji(7zH-^ — A^. In the large L limit we take jl — SttHGAl. The dihedral angles <i>„m = $ 
of a regular tetrahedron are given by cos($) = —1/4. Therefore H84|l becomes 



'J'qW = C5 exp <( -— 22 "(«m)(p<3) iJnm " JL){jpq - j l) + ^^ E ■^"™ ( ' (^"^^^ 

{7im){pq) [n.in) 

To respect the symmetry of the sphere, the covariance matrix a(„m)(pg) of the gaussian can depend 
only on three numbers 

Ot{nm){pq) = "l a{,im){pq) + "2 5(^nm)(pq) + "3 &(nm)(pg) (103) 

where 5(^nm){pq) = 1 if (nm) = {pq), ai^nm)(pq) = 1 if just two indices are the same, and 6(„m)(pg) = 1 if 
all four indices are different, and in all other cases these quantities vanish. We will use this notation, 
namely a(i2)(i3) = ai, a(i2)(i2) = "2, a(i2)(34) = "a repeatedly. 
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The component of the state H84|l that matters at first order in A is thus completely determined 
up to the three numbers ai,Q;2,a3, and the constant C5. This amounts to select a vacuum state 
which is a coherent state peaked both on the background values of the spins (the intrinsic geometry 
of the boundary surface), and on the background values of the angles (the intrinsic geometry of the 
boundary surface). See 10 for a similar construction in 3d. 

For clarity, let us stress that we are not assuming that the boundary state has components only 
on the five-valent graph considered. What we are saying is that only this component of the boundary 
state enters the expansion to first order in A that we are considering. This has generated a certain 
confusion in public discussions of preliminary versions of this work. 

5.1 First order graviton propagator 

We have now all the elements needed to compute the expression (|5^ . Inserting IjlQll) , (|102|l and l(7H| 
in (|52(l we obtain a completely well-defined expression for the propagator. As a first step towards 
the analysis of the resulting expression, we choose the points x and y to be two distinct nodes of the 
boundary spinnetwork. Equivalently, these can be thought as points located, say, in the centers of 
the corresponding dual tetrahedra: in the theory, of course, position is not determined with better 
precision that the individual "atoms of spaces" described by the individual tetrahedra. We consider 
the ten by ten matrix G(L) formed by the "diagonal" components of the propagator 

G(i)(,,.)(,,) ^ Gf^''ix,y) <-'nf > n<-'n<-', (104) 

where n^'^' is the normal to the triangle ty. Since all ten triangles have the same background area, 
|n('^'| = \n\ = S-KhGJL for large jl, we can write G(L) = G(L)/|n|4 = G(L)/(87r?iGJL)l By 
symmetry 

G(-Zj)(y)(/cO = Gi(L) a(y)(j,;) + G2{L) 5(ij)(ki) + Gz{L) b(^^j^ki)- (105) 

This expression depends of course on L, which determines the distance between the points considered 
and the angles between the directions considered. 

Before computing this quantity in the background independent theory, let us compute it in con- 
ventional linearized quantum general relativity, for later comparison. There are two ways of making 
the comparison. One is to compare the propagator obtained in the background independent calcu- 
lation with the linearized-theory propagator G^^^^j.^^^^{x , y) , where x and y are in the center of the 
corresponding tetrahedron. The other is to compare it with the quantity obtained by integrating 
^iinoarizcd('''' v) ^^ ^ ^^'^ 2/' ^^'^'^ ^^^ entire tetrahedron that they represent. The difference is a numer- 
ical factor that is not relevant for us, and we choose here the first option. In a flat background metric, 
two points in the center of adjacent tetrahedra, in a surface with the boundary geometry chosen, are 
at a distance \x — y\q — L/A. If the four indices i,j, k, I are all distinct, it is easy to see that n**^' and 
^(ki) j^j^g orthogonal; then the propagator is easily computed to be 

On the other hand, the components Gl™s?'"?'''* and G'/"??[m°'^ are vacuum expectation values at fixed 
"time" : the first is the fluctuation of the area square of a triangle, and the second is the vacuum corre- 
lation between the fluctuations of the area squares of two adjacent triangles in the same tetrahedron. 
These are also proportional to i"^. We can therefore write 

glincarizcd,^N ^ ^J;^ y^ ^^qj^ 
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where W is a numerical matrix, with the same symmetry structure as in H1Q5|I . For instance, 
yV(i2)(34) — i &s in H1Q6(I . while the others projections are easily obtained from the linear theory. 

We now compute the matrix G(L) in the full theory. Since this is a diagonal term in the propagator, 
we can use (f7H|l and (1^ reads 

GiLh,)iki) - ^^^E^W ((8^?iG)2j,,fe- + 1) - \n\'){{8nnGfjki{jki + 1) - \n\^)^M (108) 

The terms jnp come from the background 5°"^ and are equal to the square of the area of the face, 
namely to {SnhGJL)'^, for large j. Inserting (|126|l and H1Q2|I we have, to first order in A 



= ^-tE n dim(j„™) (x 

' jnm \n<7n / 

Bijnrn) C5 CXp < ^ ^—a(^„rn){pq)Unrn " JL)iJpq - j l) + «$ ^ ■?'"" f 
I ^^^ n<rn J 



G(i)(y)(fe;) = ^— zJ I n dim(j„™) I Oy (jy + 1) - ji) ijklijkl + 1) - ii) 

(109) 



where we have used the Einstein convention in the exponent. Since we have assumed that j^ is large 
and the vacuum exponential peaks the sum around jl, we can discard the +1 in the parenthesis. We 
expand the summand in the fluctuations Sjij = {jij — jl), and keep only the lowest term, assuming 
that the gaussian suppress the higher terms. This gives 

' jnm \n<m / 

B{jnm) C5 exp < -77—<^(nm)(pq) Sjnm Spq + i* E ■?""» f ' 



(110) 



We assume that the dimj terms vary slowly over the range where the gaussian is peaked, and can be 
considered constant. Let us absorb C5 and these constants in a factor A/s . 

G(i)(y)(fci) = A/'s— E '^•?"™ ^P1 ^iJnm) exp i --r— a(„m)(pg) ^jnrn Spq + i$ ^ jnm \ ■ (HI) 

■^^ hm [ ^^ n<m J 

We change summation variable from the spins to the fluctuation of the spins 

G(i)(,,)(fe,) = AA54 E ^^^^ ^J'^' ^(^'i + '^^""O e-^"<'"")<-''-'''""'^-+'* S""^ '- . (112) 

The sum can be approximated with a gaussian integral. The rapidly oscillating term expi$ "^nm Jrim 
tends to suppress the sum. To evaluate it, we need the explicit form of B{iL + Sjnm) in the large 
j regime, discussed above. Since the sum (|112|l is peaked around jnm = Jl, let us expand the lOj 
symbol around this point. To second order around jnm = Jl, the Regge action reads 

^'Rcggc\Jn7n) / ^ Jrun i 'Z^ {nnn){pq)^Jmn^Jpq, V-^-^*^/ 



where, introducing the "discrete derivative" q''' = f{j + 1/2) — f{j), we have defined 

d^mnijrs) 



G(mn)(pq) q- 
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(114) 



Thus, around jnm — JL, l|85|l gives 

(115) 
where tr is the regular four simplex (for which k^ = 1), which is the only non-degenerate four-simplex 
with these areas [HI]- The key observation is now the fact that the rapidly oscillating exp{i<i> J2nm Jnm} 
term in the second term of this expression cancels with the rapidly oscillating term in H112(l . Therefore 
the second term of the last expression contributes in a non- negligible way to the sum pi2|) . The first 
term is suppressed (by the rapidly oscillating factor exp2i$^^^ j„m) and it is reasonable to expect 
that so is the degenerate term D{jnm) because this corresponds to 4-simplices with different angles, 
and should be dominated by different frequencies. This is the cancellation of the phases that we 
mentioned in Section 2. Therefore (|112|l becomes, keeping only the first term 

G(L)(„-)(,;) =AA^4 E ^^^^ ^^fc' e-^^'-)<-)'^-^^-e-*"<""'<-'''^"'"'^-, (116) 

^^ Sum 

where we have absorbed some constant factors in M^. This factor, which contain the constant C5 of 
the state, is determined by the WdW condition ((SSI)- This is simply given by the same expression 
without the field operator insertions, namely 

l^TVj^y^ g-5G(„™)(p,)5j„™<5ip,g-27r"("")(P9)'5j"™'5>9, (117) 

Thus using the form (|54ll for the two-point function, we have 



G(i)(ij)(A;() = -J ^ •^"'" "-]—— 1 T- ■ (118) 



This expression for the 2-point function has the same structure as equation (|ll|l . Approximating the 
sum by a gaussian integral gives 

G{L) = \{rL^a + iGy\ (119) 

Jl 

We only need to evaluate the derivatives H114(l of the angles with respect to the spins. 4d geometry 
gives (see the Appendix A) 

^ ^.r^d^nm STThG (7 \_87ThG 

G(nm){pq) = GT^nG -^ = I -a(„m)(pg) - 'dd(^nm)(pq) " 4tl(„,„)(pg) I = ^^ l^(nm)(pq)- (,12Uj 

The ten by ten matrix /C has purely numerical entries. From the relation between areas and spins, we 
have JL = VSL^ / {3271110) . The jl factor that combines with the one in front of a in H119|l to give 
the crucial overall 1/L^ dependence of the propagator. Finally (|119|l reads 

G(i) = -|^(a-KzV3/4/c)"\ (121) 

This is the value (|106I107|I of the propagator computed from the linearized theory, with the correct 
l/\x — yp spacetime dependence. The three numerical coefficients of the matrix a are completely 
determined by a = Att'^/VS W"^ - iV3/4 /C. 
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Figure 1: The boundary graph T and the boundary spin network s. 

6 Second order 

Next, we consider the second order term in the expansion in A. We focus on one of the contributions, 
leaving the general case for further developments. The case we consider is determined by a boundary 
graph Fg illustrated in FigQ It consists of two tetrahedral spin networks connected by four links. 
Denote the nodes of the first tetrahedra u„ and the nodes of the second one dn, n = 1, ...,4 (u for 
"up", d for "down"). The links of F are the six links ZJJ,„ with n < m oi the "up" tetrahedron, the six 
links /,^,„ of the "down" tetrahedron, and the four side links /^, connecting w„ with dn- 

In order to illustrate the geometrical meaning of this boundary graph, consider the four-dimensional 
triangulation A illustrated in the first panel of Fig.|2J formed by two four-simplices (say the "up" one 
v"P and the "down" one v'*°™'^), joint by a tetrahedron T. If we interpret the vertical axis as a "time" 
axis, the triangulation A represents the world-history of a point d opening up to a tetrahedron T and 
then recoUapsing to a point u. This is the 4d analog of the 3d and 2d cases illustrated in the other two 
panels of the figure. In the 3d case, we have a point opening up to a triangle and then recoUapsing; in 
the 2d case, we have a point opening up to a segment and then recoUapsing. Label the four vertices of 



i-d 



3-d 



2-d 



t=T 



t = 



t = -T 




Figure 2: The spacetime triangulation A together with its 3d and 2d analogs. 
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t = 



Figure 3: The spacetime triangulation. 

T with an index n = 1, ..., 4, caU V^^^ (respectively Inm) the triangles formed by the vertices of T n, to 
and u (respectively n, m and d), and In the triangle opposite to the vertex n of T. The four triangles 
II bound r (Fig. EJ. 

Now, consider the boundary of the triangulation A. The eight boundary tetrahedra and the 
tetrahedron T are represented in Fig^ This is formed by four "up" tetrahedra u„ (the ones bounding 
t;"P, except for T) and four "down" tetrahedra (i„ (the ones bounding ti^own^ except for T). The two 
upper tetrahedra u„ and Um are separated by the triangle t"™- Similarly, the two lower tetrahedra 
are separated by a triangle i„m. Finally, each upper tetrahedron is separated by a lower tetrahedron 
by a triangle t„. 

The dual of the boundary of A is precisely F, defined above, and illustrated in FigQ the four 
upper nodes of F correspond to the four upper tetrahedra; the four lower nodes of F correspond to 
the four lower tetrahedra. The six links joining the upper (lower) nodes correspond to the six upper 
(lower) vertical triangles i"„ {inm)'^ the four vertical links correspond to the four triangles i„ bounding 
T. In conclusion, F is the dual triangulation of a 3d surface that can be viewed as the boundary of a 
spacetime region formed by a point expanding to a tetrahedron and then recoUapsing to a point. 

Let us denote j"m, j^™, js„ the twenty simple 5*0(4) irreducible representations associated to the 
20 links Inrm^nrm^n of F, and iniin eight intertwiners associated to the eight nodes u„ and d„. The 
set s — (F, j,",„, j^,„, j^, i", ij^) is the boundary spin network we consider in this section. 

The boundary function fs{4') determined by this spin network is 

is{c^) = E n ^^,:'" ^S:"" (122) 

This is a monomial of order eight in the field, and is an observable in the group field theory. Its 
expectation value is given by H59(l . At order A^, this gives 

^W^2(^/°'^^^('^K/'^') ''^'' ^^^^^ 

The Wick expansion of this integral gives two vertices, say w"p and v'^°'*'^ and nine propagators. 
Namely the Feynman graph 
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Figure 4: Decomposition 



^ 



Let us focus on the particular term of the Wick expansion obtained contracting the two vertices with 
one propagator, and contracting the other four legs of w"p (respectively w*^""") with the four "up" 
nodes m„ (respectively, with the four "down" nodes d„) of s, as represented in FigEl This is the term 



W'[s] 



A2 



2(5 



02 









n=1.4 






(124) 
where repeated representation indices are summed over. This expression still contains many terms due 
to the summation over the permutations in the propagators. To find the dominant contribution, recall 
that each closed sequence of deltas in this expression is interpreted as face. Each face contributes with 
a factor equal to the dimension of the representation. The dominant term for large representations is 
therefore the one with the largest number of surfaces. A short reflection will convince the reader that 
this is the term in which the faces is the spinfoam which is dual to the triangulation described above. 
See Figure 3. That is, the dominant term is 



W[s 



A2 



2(5 



W 



n(*«KBc)(i^|JBC>dim(j„) n '^™(-^"™)^™(-^«™) l3iJnrn,Jn),B{j^^,jn)- 



\n<7n 

If we chose the mtertwiners as above, this reduces to 



W[s] = 



2(5!) 



]^dim(j„) Jl dim(j;;„)dim(j,t„0 S(jm„, jn), S(j;tm, jn)- 



(125) 
(126) 



This is the dominant term of the connected component of the amplitude for the boundary spin 
network considered, in the limit of large representations. This is the expression we will use within 
equation (|5^ . 
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t= T 



t = 




Dual spin network 




Figure 5: On the left, the triangulated spacetime. On, the right, its dual drawn around the Feynman 
graph. 



The boundary spin network is determined by the quantum numbers (i^m,i^,*"'"^)- The spins j"j^ 
and j^ are the quantum numbers of the areas of the triangles t"^ and t^ respectively. The intertwiners 



,-u,d 



are the quantum numbers of the dihedral angles between the triangles i„ and i, 



u,d 



6.1 Vacuum boundary state 

We now have to extend the choice of the vacuum boundary state \E'q[s] considered in the previous 
section to the larger spin networks considered here. For simplicity, let us again fix the intertwiners to 
be Barrett-Crane intertwiners, as above. We need therefore to select the function 



*qW 






(127) 



We want this function to be peaked on background values {jnm ,jn ) of its arguments representing 
a given background geometry q. Let this background geometry be the boundary triangulation by it 
being the boundary of the region of R^ formed by the two pyramidal four-simplices having for basis 
a regular tetrahedron T with side of length L', and height T. 

For this, let us consider a Euclidean 4d space with cartesian coordinates t, x, y, z. Consider a 
regular tetrahedron T in the t = surface, having side length L and center on the origin of E. 
Consider the two points u = (T, 0,0,0) and d = (—T, 0,0,0). The two 4-simplices bounded by u 
(respect d) and T, taken together, define a compact 4d region. Call S the boundary of this region. 
E is a triangulated 3d metric surface. Let q be the geometry of S. This is the classical metric on 
which we want the state ^q[s] to be peaked. The barycenters of the eight tetrahedra sit on a sphere 

of radius R' = \ (^Y + (ItTt^O^ ^'^d the barycenters of the two tetrahedra sharing the same face 

of T are at a distance D' = -j from each other (see Appendix A). These two barycenters are the point 
X and y where the 2-pont function is computed. In order to match the second order calculation with 
the first order one, we want the geometrical relations between these two points to be the same as in 
the first order calculation. That is, we want R' and D' to be the same as in the first order calculation. 
For this, we have to choose (see Appendix) 



2' 



L' 



I- 



(128) 
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Figure 6: Central tetrahedron (a). Lateral tetrahedron (5). 



However, instead of fixing L' and T in this manner, let us keep them independent for the moment. 
Intuitively, L' represents the "spatial" extension and 2T the "temporal" duration of the spacetime 
region considered. 

The area of the triangles of the triangulation is easily computed from elementary geometry. The 
area of each of the triangles t%_ is 



V3w2 



Ar = —L 



while the area of the each of the triangles tj^^ and tj^^ is 



•^'-4v'^+^^+f^'^' 



(129) 



(130) 



Equation (|71|l gives us immediately the background values of the spins of s. Since we are interested 
in the large j regime, we have 



s(0) _ 



V3 2_ . 

L = ]L, 



•u(0) ^ -diO) 
Jnrn Jnrn 



32T:hG 
1 



8,.Gl^^V^ + ^' + ^^'^ 



]TL- 



(131) 
(132) 



We do not give here the explicit value of the background dihedral angles ^j"'. This value can be 
obtained by elementary geometry and plays no role in the following. 

The spin network s'"' — {Tg, jnm , jnm , jn ) represents the discretization of the geometry q on 
the triangulation A of the metric surface E (where we are disregarding the intertwiners). It will be 
our reference background spin network. As above, we choose a Gaussian peaked on s^*^' . Writing all 
spins in a single vector ji = ( j "„ , J^™ , j^ ) , we have 



where we have used aw 



and the Einstein convention. 



(133) 



V^fW 



6.2 Second order graviton propagator 

Inserting 1126(1 . ((133() and 1(78(1 in ((52(1 we obtain a well-defined expression for the propagator at 
second order in A. Choose the points x and y to be nodes ui and rfi of the boundary spin network. 
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Equivalently, these can be thought as coordinates located in the corresponding dual tetrahedra ui and 
di. We thus consider the projection 

G(L', T) = G^^-'^(x, y) UaUb n^n^, (134) 

of the propagator. Notice that \n\ = |n| and that Ua and fia are orthogonal in the background metric. 
We write G{L',T) = G(L',T)/|n|4 as above. 

Let's evaluate this quantity in linearized quantum gravity. In the background metric considered, 
the two points have the same "spatial" coordinates and a "temporal" separation T/2 (the distance 
between the center of the two tetrahedra) . Hence 

Qlincarizod/^/^yN ^ (135) 

which is the same as (|107|) if (|128ll holds. 

We now compute this quantity in the full theory. Using H78(l . H52|l reads 

GiL',T) = ^-±-^Y.W[s] ((87r;iG)2jT2(jT2 + 1) - H'mTrhGr jUJi3 + 1) - I^H^qW- (136) 

Inserting (|125|l and (|133|l we have 

G(i',T) = -LJl^J^U^Mji) (jT2(jT2 + 1) - ijTLf) (jfaOfa + 1) - Utl?) 

S(j;™, fn) BU^m.fn) Cs exp {-an- {ji " jDUi' " J^ + ^'^fji } ■ (137) 
Steps analogous to the first order case lead to 

^lt)-j Jtl sj, 

where a certain number of constants, including Cg have been absorbed into As- 

We now use the asymptotic expression (|85|l for each B, obtaining a total action given by the sum 
of Regge actions for the two 4-simplicies, S"^ + Sf^^ . In this resulting action we can treat all 
spins as independent variables, and expand it around jtl and jl 

SRegge{j^,n,jin,fn) = ^^nljnm + ^i^^Jn + &£m + 4>l? jl + \GiV SjlSj, , (139) 

where 0„ and (j)nrn, are the dihedral angles of flat 4-simplices with the given boundary geometry; and 
the second order "discrete derivatives" are 

Gw = (^-^S^) ■ (140) 

This matrix can be computed from elementary geometry, as we did in Section 5.1 (see also the 
Appendix). Being a derivative of an angle with respect of an area, for dimensional reasons, Gw 

should scale as the inverse of yj^ j\, . It is therefore convenient to define the scaled quantity 

^"' (141) 



jf'jP 
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We assume, as before, that the only surviving term in the sum l|188|l is the one in which the 
exponential in the cosine matches the phase of the boundary state. Notice that this happens because 
the linear terms in the expansion of the Regge action sum up, giving the dihedral angle of the boundary 
of the 4d region, which is precisely the sum of the dihedral angles of the two 4-simplices at the faces 
of T. That is, ^S^j = $1°„\, but 2^1°^ = $1°\ Thus, we can rewrite (fli^ as 



G(L', T) = fXf'^ V <5j- <5jf3 P2 e-<Sn...^+sL..^+'^^i)e-^u'in-,l°')U,-,l?')+^'^^°^n (] 

2(5-) JTL Sj: 

The first order term of the expansion of the Regge action cancels the phases in the state, leaving 

2i5irfTL 



G{L',T) = '^^gf " E ^^12 ^J'^ -'""' ^^--^-°'»^-^"'^-^^"' '- ^^" . (143) 



The normalization factor is determined by H50|) . Let us first provisionally assume that Cg is the only 
non- vanishing Cr coefficient in (|54|l . Then 

•^8^^ = ^ ^"i nV ' E e""""'^"^^''' e-i'^^'' *J' ^^'' . (144) 

We introduce the matrix 



Aw = 25,z- + tGw - V jf^J-r (2a + *r)„, = V^j"^jj"^ ^/f (145) 

and the vector Sj = {dji) = {Sj^,Sj^, Sj^). So 

^(^' ^) = '^ff)t-f' E ^^i'^ ^^13 e-^-^^-^^^ (146) 

^8-^-^^^Ee-^^^^-^. (147) 

And approximating the sum with gaussian integrals gives 

GiU,T)^^i^-!^^^^{A)-!^. (148) 

Under the condition H128() . G{L) = G{L'{L),T{L)) is proportional to 1/jl, as in the first order 
calculation Thus we recover the expected -p- behavior of the linearized theory. The result is not 
proportional to A. 

What if Cs is not the only non-vanishing Cr constant? To illustrate this more general case, let 
us assume that only C5 and Cg, namely the two terms considered at order one and two are non 
vanishing. Then, calling G5(L) and Gs{L) the 2-point function computed with the sole C5 and Cg 
terms, respectively, we have that the 2-point function has the structure 

G(i) = AC5a5G5(i) + X'CsasGsiL), (149) 

for appropriate constants 05 and «§; while the WdW condition (|53|l reads 

1 = ACsas + A^Cgag. (150) 
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Posing Gs{L) = G5(i) + As{L), and k = {Cg,af,) / [C^a^) we can write, to order A 



G{L) = 



G5(L) + AfcG8(L) 



1 + Afc 
= G5(L) + AfcA8(L). (151) 

Therefore, if the second order term gives the same result as the first order term, the overall normal- 
ization makes it irrelevant. If, on the other hand, it gives a different result, then this appears as a 
correction of order A to the 2-point function. 



7 Conclusion 

We have computed a first and a second order term in the expansion in A of the diagonal components 
of the graviton propagator, in a large distance regime, starting from a background-independent for- 
mulation of quantum gravity. The result has the l/\x — yp dependence on the distance expected from 
the linearized quantum theory, the expected dependence on the physical constants, and the numerical 
proportionality constants can be fixed as a condition on the semiclassical boundary state. The main 
tool we have used is the definition of general covariant n-point function, given in 149|) . 

Many issues remain open. Among these are: the calculation of non-diagonal terms in the prop- 
agator 021; the precise physical interpretation of the two expansion parameters, A and 1/ Jl |12]; 
the physical interpretation of the numerous subdominant terms, and their relation with the rela- 
tivistic and the quantum corrections to the Newton law ^; a clarification of the implementation of 
the SO{4) -^ SO{3) gauge fixing at the boundary @H1 and of the overall Lorentz invariance of the 
formalism |49| : a precise formulation of the physical meaning of the observability of the boundary 
geometry (which is a partial observable |12ll5()llCTl52j ): the possibility of computing general covari- 
ant n-point function for models without expansion parameter A; the full exploration of the dynamical 
WdW condition (|53|l on the boundary state. These issues be discussed elsewhere. 

In our opinion, the interest of the calculation is not so much in the final agreement with the 
linearized expression or in the details of the model used, but rather in the fact that it shows how 
some low-energy quantities with a transparent physical meaning can be computed, starting from the 
abstract context of a background independent formalism. The specific choices of ingredients used 
for the calculation is therefore of lesser interest, in our opinion, that the display of the feasibility of 
calculations of this sort. 



We thank Davide Mamone for some important suggestions, and Dan Christensen for his interest 
and some important inputs based on the numerical analysis of the problem that he is developing. We 
also thank Emanuele Alesci, Federico Mattei and Massimo Testa for several discussions on this topic. 
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Appendix 

A Regular simplices 

A.l Elementary 4d geometry 

We collect here some simple geometrical formulas used in the text. An equilateral triangle of side 
L has area A = ^L^- An equilateral tetrahedron of side L has volume V3 = -rh^L"^ and height 



h = \ iL. The barycenter of the tetrahedron is at a distance d — j — :t^L from a face. 

A regular 4-simplex of side L has 4- volume V4 = g§L'^- The dihedral angles & of the 4-simplex, 
defined as the angles between the outward normals to the tetrahedra, satisfy cos 8 = —1/4. The 
center of two tetrahedra are at a distance D = j from one another and at a distance R = I— from 
the center of the 4-simplex. 

Consider an equilateral tetrahedron T in the t — plane, with center at the origin and side L' . 
Consider the two points u and d with coordinates t — T and, respectively, t — —T, and spatial 
coordinates in the origin. Fix one of the faces of T, and consider the tetrahedron defined by this face 
and the point u. The height of this tetrahedron joins the point u with the center of the face, which 
is on i = and at a distance d = —^L' from the origin. The barycenter of a tetrahedron cuts this 

height in the proportion 3:1, hence has coordinates t = -j and \x\ — j—^L'. It is therefore at a 

distance R' ~ \ K^Y + (ItTk^O^ from the origin and at a distance D' = ^ from the barycenter of 
the tetrahedron defined by the face and the point d. 

If we demand that D' = D and R' = R, we obtain T = ^ and L' = \/^L. 

A. 2 The variation of the dihedral-angle with the area in a regular 4d sim- 
plex 

Consider a 4-simplex in Euclidean 4d space. Label its vertices with an index i — 1, ...,5. Let ti be 
the tetrahedron opposite to the vertex i and fii a normalized vector normal to this tetrahedron and 
pointing outside the 4-simplex. Then the (external) dihedral angle $ij is given by 

cos$ij = fii ■ hi. (152) 

We can fix the geometry of the tetrahedron by giving the ten lengths Lij of the ten sides {i,j). 
Alternatively, we can fix the geometry of the tetrahedron by giving the ten areas Aij , where Aij is the 
area of the triangle opposite to the two vertices i and j. If we do so, the dihedral angles are functions 
^ij{Aki) of the areas. We now consider the ten by ten jacobian matrix 

J<{ij)(lm){Akl) = -TTT—- 

In particular, using the notation of (|120|l . let 

1 .. a$, 



L^^'^^'^^'^'^ - Q^^'^ 
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(153) 



be the value of this jacobian on a regular tetrahedron of equal sides Ly = L and therefore areas 
Aij = Al = ^L^ ■ Here we present a simple way of computing it, which can be easily generalized to 
more complicated situations. 

We start from the relation 

sin$yT,T, = ^FAy, (154) 

where V is the 4- volume of the simplex and T,; is the volume of the tetrahedron ti. This relation 
can be derived as follows. Pick one vertex, say « = 5 as the origin of the coordinates, and choose 
coordinates such that the four sides originating from this vertex define the coordinate axes and the 
unit length. Let gab be the Euclidean metric in these coordinates. The triangle (3, 4), bounded by the 
vertices (5,1,2) is characterized by the bivector w^^ = 1/2(5^(52 — (^2 '^i ) and its area is 



^34 = |W12| = 1/2^/511522 - (.912)^- 

The normal to the tetrahedron ta is given by n° = 1/6(5!^ and the volume of the tetrahedron is its 
length, namely T^ — |n°p = 1/6^ 5"°. Equation H152|l gives 

cos ^ab = I II I = A-^ . 155 

Hence 

(sincI>34T3T4)2 = (1 - COS^ cI>34)(T3r4)2 = 6^(533^44 „ (^34)2) ^ 6^ det <?(5ll522 - (.912)'). 

which is H154(l . because the 4- volume is 1/24 ^dct g. (The overall 3/4 factor can be easily checked by 
taking a rectangular 4-simplex of side 1.) 

Varying 1)154(1 . we have 

^ ,^ '^VAi, fdV dAi, dTi dTA . ^ fdV dA,, dT, dTA ,_^, 
cos*.,*, . j^ (_ + -j^ - _ + ^j . „„*, (_ + -j^ - _ + ^) ^ (156) 

We now consider this relation for a regular 4-simplex. In this case, cos$y = —1/4, therefore 

^fdV dAi, dT, dTA , , 

First, let us vary only one single length L^j, keeping the other nine lengths fixed. If we vary all 
lengths, then by dimensional analysis 

dV dL 

If, instead, we vary a single side, we have by symmetry 

dV^j_dL _^_._ __„^ 
V - 10 L' ' ~-^' ^ - ^ r ■ il58j 

because there are ten sides in a 4 simplex, three in a triangle and six in a tetrahedron. If we now 
vary the length of one single side in the 4 simplex, then the variation of the angle <I>.y will be different 
according to the relative position of the angle and the side, that is, according to whether or not the 
side belongs or does not belong to the triangle Aij or the tetrahedra Tj. Explicitly, 1)157(1 and ((158(1 
give 

d'^^, = -VT5 f 4 + l^v - I(e^ + ^j)) ^- (159) 



dAij _^dL 
Aij L 


dTi 
T, 


dL 
-^ L- 


y symmetry 






dAij 2 dL 
Ay 3 L ' 


dT, 
T, 


3dL 
~6T 



10 3 " 6' ■" J L 
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where e^- = 1 if the side varied is in the triangle Aij and zero otherwise and where e^ = 1 if the side 
varied is in the tetrahedron Ti and zero otherwise. 

With respect to a given angle <i>y, there are three kinds of sides. One side (i,j), six sides {i,k) 
with k y^ j, and three sides (fc, I) with k, I ^ i,j. In the first case, the side does not belong neither to 
Aij nor to Ti or Tj . Hence 

In the second case, the sides belongs to one tetrahedron, but not to the triangle, hence 

In the third case, the side is in the triangle and both tetrahedra 

In these equation and below, different indices are assumed to take different values and there is no sum 
over repeated indices. Summarizing 

^^_2/^i ^£!i^i/31 ^^_i/31 HfiS^ 

dL,j MbL' dL,u 2V 5L' dLui i^l b L' ^ ' 

The next step is to compute the variation of the sides when one area varies, at constant value 
of the other areas. By symmetry, if only the area of one triangle changes, then the variation of the 
three sides belonging to this triangle will be equal to themselves. Let (Jitriangie be such variation. 
Similarly, the variation of the six sides that have only one vertex in common be triangle will be equal 
to themselves. Let SLside be such variation. Finally, let i5I/opposite be the variation of the single side 
opposite to the triangle. 

In an equilateral triangle the three variation of the the sides contribute equally to the variation 
of the area. Consider a triangle adjacent to the triangle that varies. Its variation is proportional to 
(^■C'triangie + 2(5Lsido and vauishcs. Hence (5isidc = ^ l/2(5£triangic- The same argument applied to the 
triangles with no side in common to the varying triangle gives ^iopposito = — 2(5Lside = ^^triangie- 
Therefore we need only to compute (JLtriangie- This is easy, because for an equilateral triangle 5 A/ A = 
25L/L. Hence 5L = LjliAjA = ^bAjL'^. Summarizing: 

dU, „ 1 1 dLi, 1 1 1 dU, „ 1 1 ,_^, 

jL — 2 — — = — — = 2 — fl64) 

dAki V3 L' dA,k 2 V3 i' dA,j ^ L ^ ' 

The last step of the calculation is to combine (|153|l and Hlt)6|l 



a$y _ Y^ a$y dLr, 



dAki ^—' dLnni dAki 



(165) 



where care should be taken, for each term of the sum to use the appropriate value depending on 
whether or not there are repeated indices. This is straightforward and gives 

(166) 



9$,j 


1 1 


d<^,j 


7 1 


5$,y 


^ 


dAki 


V5i2' 


dMk 


2 V5£2' 


dA,, 


X/5L2 
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In conclusion, we can write H153|l as 

;^/C(,,)(,,) = -i^ ga^ - 4^^ - 96^) (167) 

where 5f- = 1 if the two couples of indices are the same, aj^- = 1 if two and only two indices are the 
same, and &^j = 1 if all four indices are different, and in all other cases these quantities vanish. 

B Boundary intertwiners 

At the end of Section IT^ we have mentioned the difference between the space of the simple 50(4) 
intertwiners and the space of the SU{2) intertwiners. The difference shows up in the linear structure 
of the states. Suppose we decide to pair the four faces of the tetrahedron differently, and to represent 
the intertwiner in terms of the virtual link k = ju" for a different pairing. In both cases, the linear 
properties of the space of the intertwiner allows to express a virtual link as a linear combination of 
virtual links of a different pairing, but the linear structure is different. In fact, in the SU{2) case, the 
recouping theorem gives (see for instance (A. 65) of J21j also for the notation.) 



J2. ; , J3 



>^ ^?{jM:j} 



Jl J4 




(168) 



On the other hand, in H73() the simple intertwiner labelled by j is in fact formed by a couple of 
intertwiners, one for the left and one for the right component of 50(4), having the same spin. Namely 



J3 J2 




(169) 



where the first diagram represents the self-dual and the second diagram the anti-selfdual components 
of the representation. Therefore in the 4d case a simple link in one pairing is equal to a sum including 
non simple links (when k =/= I) links in another pairing. That is 

k,l k 

The reason of the discrepancy between the linear structures in 1] 169(1 and (|170|l is not entirely clear 
too us. This discrepancy, on the other hand, does not affect the computations in the theory or the 
interpretation of the boundary states in the model we are considering. The reason is that the only 
intertwiner appearing in the spinfoam sum is the Barrett-Crane intertwiner, which decomposes into 
simple virtual links for any pairing: 



Jbc 



= E (2^-+i) y< y< = E (2^+1) A A (1^1) 
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Therefore if we identify the intertwiner ij , with the LQG intertwiner jj , , we obtain simply and 
consistently equation lf75|l. 



C Non-diagonal terms of the graviton operator 



The non-diagonal terms can be computed from 

E'\n)Ef{n)\s)^hG\s,,ij) 



(172) 



where Snij is the spin network obtained by adding an infinitesimal link of color 1, joining ("grasping") 
the links / and J infinitesimally close to the node n. This is a spin network with the same graph and 
the same spins associated to the spins, but a different intertwiner at the node. Assume for simplicity 
that the pairing chosen to label the intertwiner is such that the two grasped edges are paired together. 
Simple recouping algebra gives 



Tet 



E'^n)E^ {n)\j) ^ hG 



J2 J2 1 



\j) 



^(jl,j2,j) 

where and jiand J2 are the spins of the links / and J. This gives (T2;, pg 258) 

E'\n)Ef{n)\j) - HG ^ b^J + 1) - Ji(ji + 1) - J2(j2 + 1)] |j). 



(173) 



(174) 



In the case in which the grasped links are not the paired ones, the action of the operator can 
simply be obtained by changing basis with the recouping theorem. This gives 



E'\n)E^{n){x)\j)^hGj2c';\k) 



where 



c^ = nGY: 



ii h I 

J3 ji 3 



Tet 



J2 J2 1 



ji h k 



&iJl,J2,j) 

Here j'l and J2 are the spins of the two grasped links and j^ and j^ the other two 
In conclusion 

{k\E'\n)E^in)\j) = hGj{j + l)S,k 
when I = J; here j is the spin of the grasped link. 



Tet 



{k\E^\n)Ef{n)\j) 



TiG- 



31 31 3 

32 32 1 



0{iid2,j) 
for the grasping on paired links; here j is the virtual link. And 



3jfe 



\E'\n)E'l{n)\j)=hGY. 



31 32 I 
33 34 j 



Tet 



ji ii I 
h h 1 



^(ii,j2,j) 



31 h k 



(175) 



(176) 



(177) 



(178) 



(179) 



for the grasping of unpaired links. The action of E°'''{n)E^{Ti) is the immediate from E'^'^{n)E^{n) = 
n'}n'^jE^^{n)E^ (n) where rtj is the inverse matrix to n^. 

The calculation of the non-diagonal terms of the propagator, involving the action of these opera- 
tors, is in progress [T7j . 
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